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,-£h . Abstract: In this work we present the results of full Geant4 and FLUKA simulations and com- 



parison with dosimetry data of an electron LINAC of St. Maria Hospital located in Terni, Italy. 
The facility is being used primarily for radiotherapy and the goal of the present study is the de- 
tailed investigation of electron beam parameters to evaluate the possibility to use the e— LINAC 
(during time slots when it is not used for radiotherapy) to test the performance of detector sys- 
tems, in particular those designed to operate in space. The critical beam parameters are electron 
energy, profile and flux available at the surface of device to be tested. The present work aims to 
extract these parameters from dosimetry calibration data available at the e— LINAC. The electron 
energy ranges from 4 MeV to 20 MeV. The dose measurements have been performed by using an 
Advanced Markus Chamber which has a small sensitive volume. 
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1. Introduction 

The e— LINAC in use at Azienda Ospedaliera di Terni (AOT) has an electron energy range from 4 
MeV to 20 MeV which is of major interest in the field of performance testing of devices operating 
in space. The e— LINAC of AOT is intensively used for radiotherapy of the patients, however, 
with a special agreement and in cooperation with AOT we have investigated the possibility to use, 
during nights and the weekends, the e— LINAC to test particle detector devices and components. 

For this work, we set up a full geometrical description of e— LINAC and its electron beam to 
reproduce the dosimetry calibration data widely available at AOT. Using the Geant4 and FLUKA 
packages [|l]-@] allows the full simulation of the e— LINAC operations as well as to convert dosime- 
try data to parameters such as particle flux, energy spectrum and spot size at the surface of the 
devices positioned at various distances with respect to the e— LINAC head. 

2. Dose measurements 

The Advanced Markus Chamber (AMC), which is of PTW type 34045 (its main parameters can 
be found in ||5|]), was used to perform dosimetry measurements in AOT. The Advanced Markus 
Chamber is a parallel plate ionization chamber and has a small sensitive volume, 0.02 cm 3 , with 
a thin entrance window. The small sensitive volume of the chamber allows a dose distribution 
measurements in air and water, with good spatial resolution. The AMC used in the measurements 
can be seen in Figure [[]. 
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Figure 1. Advanced Markus Chamber (AMC). 



In AOT, dose measurements have been performed at all available energies, ranging from 4 
MeV to 20 MeV, with the trimmer (collimator) having a window surface of 3x3 cm 2 . The AMC 
has been placed at distances here called Source to Surface Distance (SSD). The setup arranged for 
the dose measurements is shown in Figure ||. 

The AMC dose profiles have been evaluated at various distances (z-axis) and different energies 
in the transverse direction (x-axis) by measuring the FWHM of maximum electron beam intensity 
observed at the central axis position during an exposition time of one minute. 

3. Simulation 

The Geant4 and FLUKA simulation tools are developed primarily for use in High Energy Physics 
applications but are used in many areas such as space applications, medical applications and accel- 
erator driven systems. We have developed simulation codes using the versions Geant4.9.3.p01 and 
Fluka2008.3c. 



a 
_■ 

II 

O 



Figure 2. A sketch of experimental setup used to take calibration data with Advanced Markus Chamber 
(AMC). 

The geometry used in the simulation code consists of an accelerator head and small air ioniza- 
tion chambers placed along the scanning axis. The accelerator head shown in Figure [| consists of 
internally embedded collimator leafs composed of a high density material and an applicator shield 
(trimmer) made of a high melting point Lipowitz material containing 50 % bismuth, 26.7 % lead, 
13.3 % tin and 10 % cadmium by weight. Geant4 and FLUKA use different functions for geomet- 




-2- 



rical definitions, but the description of the experimental setup and its operation is the same in both 
frameworks. 



The primary electrons in both simulation codes were generated at a certain point and oriented 
along beam axis. Afterwards they hit a scattering foil followed by the collimator leafs inside the 
accelerator head and a trimmer. In order to obtain a preferred beam at the accelerator exit as 
foreseen in its manual the distance between collimator leafs has been kept flexible for tuning. 

One of the key parameters to characterize the passage of a particle in a detector is its energy 
loss within the material. Since the particles of interest are electrons, we have used the Penelope 
physics, which is one of the recommended built-in physics converted from the PENELOPE code, 
for Geant4 and PRECIS 10, which gives more accurate results with low particle-production and 
tracking thresholds, for FLUKA. Considering the relatively small thickness for the AMC sensitive 
region, where ionization takes place, a smaller step size has been applied for that part of the detector 
region. Detailed information about the underlying physics of particle interactions with matter has 
been implemented in our simulation frameworks can be found in Geant4 Physics Reference Manual 
[@] and FLUKA Online Manual 0]. 

3.1 Beam Profiles 

The beam profile FWHM widths measurements at a given SSD were performed by placing the 
AMC at positions perpendicular to the beam axis. As the SSD increases, due to small divergence 
of the beam and mostly because of multiple scattering, wider profiles were measured. Simulated 
variations of beam profile width spreads are plotted for two mid electron beam energies supported 
by the e— LINAC, at 6 MeV and 15 MeV, in Figure |j. The shapes predicted by the two Monte 
Carlo simulations agree quite well. 

A comparison of the FWHM measured by the AMC with the corresponding simulations from 
both Geant4 and FLUKA are presented in Table I. 

3.2 Dose Profiles 

The Monte Carlo simulations directly predict the energy released into a specified volume, which 
is the AMC in this case. The relevant energy loss mechanism is the ionization in the chamber's 
sensitive volume. 




Figure 3. Geometry of simulated acceleration head. 
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Figure 4. Change of beam profile widths at various exposition distances for all particles. 



This deposited energy given in MeV in both simulations has been converted to the relative 
dose in order to compare it with the experimental data measured by the AMC. The relative dose 
values were then obtained at the central axis (x=0) by normalizing the dose values at a given SSD 
to the closest one at SSD=100 cm. 

The variations at two different beam energies, 6 MeV and 15 MeV, of the relative dose normal- 
ized to the one received at the central axis (x=0) are displayed in Figure |5| Experimental data and 
Monte Carlo simulations are in good agreement over all exposition distances SSD. A smaller en- 
ergy loss in the sensitive volume at larger distances is observed. In Table II we compare the relative 
doses measured by the AMC and the estimated ones from the two Monte Carlo simulations. 



3.3 Kinetic Energy at the Surface 

We have shown that the beam and dose profiles accumulated in one minute at various conditions 
can be reproduced correctly. Thus, we can evaluate the kinetic beam energy, the other crucial beam 
parameter, at the detector surface for a given distance. As the distance increase, the kinetic beam 
energy at the detector surface decreases while its fluctuation increases due to the energy spread. 
These two behaviors can be seen in Figure [6[ 
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Figure 5. The relative dose normalized to the central axis (x=0) at various exposition distances for 6 and 15 
MeV. 
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Figure 6. Kinetic beam energy and its fluctuations as a function of distance. 



The knowledge of the kinetic beam energy and its spread is crucially important to tune the 
correct beam parameters (flux, energy and profile) for the detector under test. 
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3.4 Conversion to the Flux 



The beam flux at a given distance can be evaluated by the following expression: 
Flux = CF x [ SurfaceHits/DetectorSurface/time ] 



(3.1) 



where, the conversion factor (CF) is the factor to normalize dose obtained by simulation to the 
experimental dose values and time is given in one minute of radiation at a given distance from the 
AMC. To study the performance of a particle detector at such an accelerator facility with high beam 
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Figure 7. Electron beam fluxes as a function of SSD. 



current may require reasonable extension in z direction to reduce the flux available on the detector 
surface to the desired level. From the fits performed in the simulations as illustrated in Figure [7|, 
one can extrapolate the expected lower particle flux at larger SSD. 
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Figure 8. Predicted electron beam fluxes as a function of SSD. 



The higher distances are necessary to reduce the beam flux for detector testing as can be seen 
in Figure |8[ For the electron beam energy of 6 MeV a proper distance to mount testing detector 
was found close to or larger than 250 cm while it is around 350 cm for the 15 MeV electron beam. 
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4. Measurement with CMOS Pixel Detectors 



In addition to aforementioned studies based on conversions from available AMC data, an experi- 
mental measurement of the electron flux has also been performed using a CMOS pixel detector to 
count single electrons. Our group has investigated the charged particle detection capabilities for 
this kind of sensors in the past years [§-11]. The detection efficiency for such kind of detector 



is 99 % for a minimum ionizing particle, so measuring the flux is a matter of counting the single 
electrons and normalizing it to the integration time and the detector surface. 

The CMOS sensor used in this case is a commercial VGA optical camera from Aptina Imaging, 
the MT9V01 1, as shown in Figure ||. It is a 640x480 pixel matrix, with 5.6x5.6 /im pixel size. 




Figure 9. Micron DAQ setup: the MT9V01 1 sensor with its DAQ board (Demo2). 

An electron passing through the sensor generates a signal shared only among few pixels (3—4 
on average), even for small pixel size as in our case. The frame rate for the CMOS detector is 30 
frames/second. The flux is computed as number of detected electrons/sensor surface/integration 
time. All the pixels belong to the 3x3 matrix centered on the pixel with the highest signal. The 
clustering uses two thresholds algorithm; first threshold, to define the seed pixel, requires the pres- 
ence of the signal 10 times the single pixel noise, second threshold, to define the neighbor pixels, 
requires 2.5 times the single pixel noise applied to topologically connected pixels. Hence it is a 
safe assumption to detect two different electrons if the related 3x3 matrices are non— intersecting. 
This leads to a measuring capability of roughly 10 % of the total number of pixels of the detector, 
which translates in our case to 3.5 x 10 5 electrons/cm 2 /s. 

The linearity of the measurement with the beam intensity has been calculated using the e— LINAC 
at electron energy of 10 MeV and at different distances. For beam intensities ranging from 25 U.A. 
up to 200 U.A., a linearity at the distances measured with CMOS sensor has been confirmed as seen 



in Figure 10 (left). The prediction of the Monte Carlo simulations compared to the measurements 
at available distances is presented in Figure [R] (right). Large error bar at 25 U.A. are caused by 
fluctuation of the beam. This value is a lower limit for the accelerator, seldom used in clinical prac- 
tice. The error bar at 140 U.A. is essentially dominated by statistical fluctuations, hence dominated 
by the Poissonian statistics where the absolute error is growing while the relative one is becoming 
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Figure 10. Predicted electron beam with kinetic energy of 10 MeV fluxes as a function of SSD. 



smaller. However, a good agreement within the errors has been achieved thereby validation the 
Monte Carlo setup. 



5. Conclusion 



The e— LINAC at (AOT) the Terni Hospital has a large experimental hall to perform detector testing 
with electrons of energies ranging from 4 MeV up to 20 MeV. This work describes full Monte Carlo 
simulation of AOT's e— LINAC system and compares the resulting characteristics of the parameters 
like beam particle flux, kinetic energy and its fluctuations with experimental measurements. The 
dose measurements were done with an Advanced Markus Chamber (AMC) at different Source to 
Surface Distance (SSD) and with one minute radiation time at each distance. To compare with the 
experimental data, the simulations have been performed using two different simulation packages. 
The simulations have shown that even with only AMC data one can extract relevant parameters 
such as beam energy, profile and flux, which are essential for detector performance studies. More- 
over an additional measurement to validate the Monte Carlo results by CMOS pixels using mea- 
surement technique of single electron counting has also been done resulting in a good agreement 
between Monte Carlo predictions and experimental measurements. During the data-simulation re- 
sults checks, a systematic study was done to select the correct particle production thresholds (ppt) 
and step sizes (ss) as well as physics options used in Geant4 simulations. The PENELOPE package 
was used (with ppt= 0.1 mm and ss=l mm) since it has resulted to provide better agreement with 
data. 
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6. Appendix 



Table 1. Beam Profile Widths [cm] as a function of SSD [cm] for All Data Set Compared to MC. 
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Table 2. Relative Dose Values as a function of SSD [cm] for All Data Set Compared to MC 
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